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T^' Abstract 

Conventional kinesin is a motor protein, which is able to walk along a 

microtubule processively. The exact mechanism of the stepping motion and 

Q ■ force generation of kinesin is still far from clear. In this paper we argue that 

* c/3 . neck linker docking is a crucial element of this mechanism, without which 

^1 the experimentally observed dwell times of the steps could not be explained 

f^ ' under a wide range of loading forces. We also show that the experimental 

data impose very strict constraints on the lengths of both the neck linker 
and its docking section, which are compatible with the known structure of 
'i> ' kinesin. 

in 

m 1 Introduction 

O I Kinesin motor proteins execute a variety of intracellular transport functions by 

(^ I transporting cellular cargo along microtubules (MTs) while hydrolyzing adeno- 

sine triphosphate (ATP) [6, 4, 7, 24]. These molecular walking machines move in 
8-nm steps toward the plus end of microtubules, turning over one ATP molecule 
per step under a range of loads [17, 23, 26]. Conventional kinesin (kinesin-1) is a 
homodimer, the monomers of which consist of a head domain (containing a con- 
served catalytic core), a stalk (through which the monomers form a dimer), and an 
approximately 13 amino acid long neck linker (connecting the head to the stalk). 
The neck linker is evolutionarily highly conserved among plus-end directed mo- 
tors [24] and appears to be crucial for motility [ ]. 

In a kinesin dimer the two stalks form a coiled coil, to which a cargo is attached 
or, during experiments, a pulling force is applied. The catalytic core of each head 
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is responsible for binding and hydrolyzing ATP [13, 14, 19]. ATP binding to a 
MT bound head has been shown to result in a section of the neck linker binding to 
the head [18, 3, 19, 22], with a consequence of positioning the remaining unbound 
section (and also the diffusing other head) closer to the forward binding site [17, 
15]. Experimental studies [ ] of this conformational change (referred to as 'neck 
linker docking') strongly support a scenario wherein an unstructured, random- 
coil-like neck linker folds onto the core as a result of nucleotide binding - paying 
a large entropic cost compensated by a large enthalpic gain (both in the order of 
50 kJ/mol or 20 ksT, where ks ~ 1.38 x 10"^^ J/K is the Boltzmann constant and 
T denotes the absolute temperature, which we set now to be 293 K). 

Neck linker docking is also confirmed by X-ray structures [12, 21], molecular 
dynamics simulations [8, 11], and fluorescence experiments [16], however, it is 
still debated whether docking is crucial to the processivity and force generation 
of kinesin, or just a byproduct of some other mechanism [10]. To address this 
question we focus on a single state of the dimer, in which one of the heads is 
bound to the MT and contains an ATP, while the other head (often called the 
'tethered head') is unbound and contains an ADP (see left cartoon in Fig. 1). To 
complete a forward step, the tethered head must find and then bind to the forward 
binding site on the MT (which is about L = 8 nm ahead of the bound head). 

In a recent optical tweezers experiment Carter and Cross [ ] demonstrated 
that kinesin can take not only forward but also backward steps. The ratio of the 
number of backward and forward steps increases as the loading force is increased 
and reaches unity at the stall force of about —7 pN (a negative force means pulling 
toward the minus end of the MT). At saturating ATP concentration and at zero 
loading force the steps almost exclusively occur in the forward direction and their 
duration (which is also called the dwell time) is around 0.01 s. Near the stall 
force, however, the steps are much slower with a dwell time of about 0.3 s. We 
note here that in any situation when two alternative processes compete with each 
other (such as the forward and backward steps of kinesin) the duration of the 
successful process will be determined by that of the faster one (meaning that the 
completion of the slower one will take a similar amount of time, but will occur 
less frequently). So at small loads the forward steps, whereas at loads exceeding 
the stall force the backward steps will determine the dwell times of the steps in 
both directions. 

The benefit of neck linker docking appears to be clear: in the one head bound 
state (described above) this conformational change positions the tethered head 
closer to the forward binding site on the MT (as demonstrated in Fig. 1). However, 
to underpin the necessity of neck linker docking in a quantitative manner, we 
examine under what conditions can the tethered head find the binding site within 
the aforementioned 0.3 s time limit when the —7 pN stall force is applied. 
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Figure 1: One head bound kinesin before (left) and after (middle) neck linker 
docking, followed by the binding of the tethers head to the microtubule (right). 
Inset shows the number of docked neck linker segments (iVd) in white and that of 
the remaining (A^^) segments of the same neck linker in black. 

2 The Model 



We suppose that the undocked sections of the neck linkers behave as random coils, 
which we simply model as freely jointed chains (FJCs). Thus, before binding 
to the MT the tethered head experiences a random walk in the potential of the 
entropic springs formed by the two connected neck linkers of kinesin. We regard 
the spatial distribution of the tethered head simply as that of the end of the two 
joined neck linkers. Each neck linker consists of A^aa amino acids (A^aa ~ 13 for 
a conventional kinesin). We use /aa = 0.38 nm for one amino acid length and 
/p = 0.44 nm for the persistence length of a polypeptide chain [20]. The Kuhn 
length is then /k = 2/p = 0.88 nm, and the number of freely jointed segments 
in a neck linker is A^ = A'^aaLa/^Kj vvhere A'aaLa is the contour length of the 
polypeptide chain [, ]. As one can see this description is quite coarse with respect 
to the number of amino acids of the neck linker: increasing the number of the 
segments A^ by 1 is equivalent to increasing the number of the amino acids A'aa 
by /K//aa ~ 2.3. 

For technical reasons we divide the neck linker of the bound head into two 
sections. The first one consists of those A^ segments that can dock to the head 
domain, while the second one contains the remaining A'u = ^ — -^d segments 
of the neck linker. The external pulling force F acts at the end point of this 
neck linker, where it joins to the neck linker of the tethered head (as illustrated 
in Fig. 1). For simplicity we neglect the width and flexibility of the coiled coil 
stalk between the two neck linkers. This omission could, however, be largely 
compensated by a slight increase in the length of the neck linkers. 



Let X, y, and z be Cartesian coordinates so that the x axis is parallel to the 
MT, and the x and y axes span the plane of the external force F. This way F has 
no z component. The angle of the force {i.e., that of the coiled coil stalk) to the 
MT depends on the details of the experimental setup, in particular, on the length 
of the stalk and the size of the bead in the optical trap. Throughout the paper we 
use a reasonable value of 45° for this angle and, thus, assume that Fy = \Fx\. In 
the experiments, usually the x component of F are reported, so by specifying that 
the stall force is —7 pN we mean F = (—7, 7, 0) pN. 

Let p%{R) denote the probability density of the end-to-end vector _R of a free 
(F = 0) neck linker with N segments. Such a distribution applies to the neck 
linker of the tethered head. As p%{R) does not depend on the external force it has 
a spherical symmetry. The formula for this distribution is derived in the Appendix. 

Let Pn{R, F) denote the probability density of the end-to-end vector R (point- 
ing from the head end of the neck linker towards the pulled end) of the MT-bound 
head's neck linker with A^ segments at an applied external force F. Note that the 
applied force breaks the spherical symmetry of the distribution. Pn{R, F) can be 
expressed with the help of p%{R) as 
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Pn{R,F)= (1) 

^N{-r ) 

where 

Zn{F) = Jp%{R)e^dR (2) 

is the partition function. 

The compound probability density pn^^n2{R^ P) for the end-to-end vector R 
of two joint neck linkers (with A^i and A'2 segments, respectively) that are pulled 
by an external force F at the joint (and held fixed at the other end point of the first 
neck linker) can then be expressed as: 



PNi,N2 



{R, F)= j PN, {R', F)p%^ {it - R')dR' . (3) 



By placing the origin of the coordinate system to the starting point of the neck 
linker of the bound head, the probability density pj^^iR, F) can be considered as 
the concentration of the tethered head at position R, given that the neck linker of 
the bound head is undocked, and if all steric constraints are neglected. 

To take the volume exclusion between the tethered head and both the MT 
(approximated as Ry < 0) and the bound head (approximated as |/^| < 2 nm) into 
account, we introduce a constraining function 

0Cm ^ [ ifi?2^ < Oor 1^1 < 2nm, 
1 1 otherwise. 



Multiplying pj^^iR, F) by Q{R), and then renormalizing it to unity results in the 
sterically constrained concentration of the tethered head: 

IPKAR,F)e{R)dR 

Note, however, that the volume exclusion is not imposed on the entire chain, only 
on its end point. Nevertheless, this steric constraint has very little effect on our 
main conclusions. 

Similarly, denoting the end-to-end vector of the N^^ segments of the docked 
section of the neck linker by L^ = {L^^, 0, 0), i.e., assuming that docking occurs in 
the X direction with a projected distance of Ld, the probability density p^^^NiR — 
Ld, F) can be considered as the concentration of the tethered head at position R 
when that the neck linker of the bound head is docked (and all steric constraints 
are neglected). After its multiplication by Q{R) and renormalization to unity, 
results in the sterically constrained concentration 

IP!,.AR-Li,F)e{R)dR 

Multiplying these concentrations at the positions R = L = (L,0, 0) and R = 
—L = {—L, 0, 0) by the binding rate constant k^ of the kinesin head, we get 
the rate at which the tethered head binds to the MT at the forward and backward 
binding sites, respectively. 

The only remaining quantity to determine is the probability of the neck linker 
being docked or undocked (disregarding the possibility of any long lived partially 
docked conformation). These probabilities must be a function of the force since 
the neck linker docks easier if it is pulled forward. Let P*(F) denote the proba- 
bility that the neck linker is docked and P{F) = 1 - P*{F) that it is not. Then 
the free energy difference AG{F) between the docked and undocked state can be 
defined through 

P*(F) -AG(F) 

— ^ = e ^bt . (7) 

P{F) 

From the measurements of Rice et al. [ ! 8] we know that if an ATP is bound 
to the head, then neck linker docking is an energetically favorable process with a 
free energy difference of AG" ~ —2k-QT. Thus, the ratio of the probabilities that 
the bound head is in the docked and the undocked conformation can be written as 

P\F) _ e '-bt llpl^[R!-Uy-.-p%{R-R!)Q{R)d.K^R 

jjiB\ ~ z — m z — z z — zr^ • ^^> 

^(^) JJp%{R')e^p%{R-R')e{R)dR'dR 
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Figure 2: The characteristic time tf needed to find the forward binding site as a 
function of the horizontal projection L^ of the docked neck linker at the —7 pN 
stall force. The four boxes correspond to neck linkers of different lengths with 
different numbers N of Kuhn segments. The various lines in the boxes correspond 
to different numbers A^d of docked segments. The rectangular area in the bottom 
left corner of each box indicates the expected range tf < 0.3 s and L^ < 3.5 nm, 
compatible with the experiments. 



The characteristic time for binding forward at force F is then: 
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U{F) 



h[P*{F)c*m F) + P{F)c{L, F)] 



(9) 



We have no reason to assume a different rate constant fcb for backward and 
forward binding, therefore, the ratio of the probabilities of forward and backward 
stepping is: 

P*{F)c*{L,F) + P{F)c{L,F) 



r{F) 



P*{F)c*{-L, F) + P{F)c{-L, F) 



(10) 
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Figure 3: The ratio r of probabilities to bind forward and backward at the —7 pN 
stall force as a function of the horizontal projection L^ of the docked neck linker. 
The four boxes correspond to different values of A^, and the different lines in the 
boxes correspond to different numbers A^d of docked segments. The preferred 
range of r > 1 and Ld < 3.5 nm can be seen in the top left corner of each box. 



3 Results 

We evaluated the integrals numerically for four integer values of A^ between 5 and 
8. A^d took values between and A^ — 1, so that we could compare the situation 
A'd = 0, where neck linker docking plays no role in the searching of the bind- 
ing site, to various neck linker docking geometries. Although for completeness 
we varied the horizontal projection L^ of the docked neck linker from 0.4Ad^K 
to A^d^K^ only values below 3.5 nm are acceptable since the kinesin head has a 
diameter of around 4 nm [ i ] . 

The results for the characteristic forward binding time tf at the —7 pN stall 
force are presented in Fig. 2. All the curves decrease monotonically on the interval 
Ld G [0, 3.5nm]. Thus, if we want to find the fastest forward step it is enough to 
look at the curves at L^ = 3.5 nm. The rectangular areas in the bottom left comers 
of the four boxes indicate the desired range of tf < 0.3 s and the allowed regime 
of Ld < 3.5 nm. 
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Figure 4: The same as Fig. 2, but at zero load, and the experimentally acceptable 
range is tf < 0.01 s. 



We used /cb = 20 s~^/iM~^ for the rate constant of MT binding [25, 9]. None 
of the curves passes through the desired area. We can make them do so by in- 
creasing fcb (which would simply shift them downwards). This can be justified 
by noticing that the k\, ^ 20 %~^ijMr^ value was measured for free kinesins in 
solution. However, in our case the tethered head is diffusing in the close vicinity 
of the MT, where electric fields are not fully screened and local attractive inter- 
actions might be significant, which can lead to an effectively elevated value of 
/cb- A threefold increase seems sufficient to make some curves just pass through 
the desired area for A^ > 6. A less feasible tenfold increase would, however, be 
necessary for A^ = 5. 

In the r vs. L^ plots (Fig. 3) the value of /cb does not play any role (see Eq. 
(10)). All the curves for A^ > 6 are orders of magnitudes below r = 1, the value, 
which would be the expected at the stall force, if each forward binding resulted in 
a forward step, and each backward binding in a backward step. Obviously, there 
must exist some rectification mechanism that makes the completion of a backward 
steps less likely {e.g., by reducing the ADP dissociation rate from the rear head) 
and increases the value of r to 1, but this is not the subject of our investigation. 
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Figure 5: The same as Fig. 3, but at zero load. 

For comparison we repeated the calculations at zero load as well (Figs. 4 and 
5). These results indicate that for many A^, N^, and L^ combinations U is below 
the 0.01 s limit imposed by the experiments, even at k\^ = 20 s~^/iM~^ Also, 
the values of r for L^ = 3.5 nm are orders of magnitude larger than 1, so with 
neck linker docking backward steps are highly improbable even without any rec- 
tification, in accordance with the experiments. Nevertheless, as r can never be 
smaller than 1 at zero load, the same rectification mechanism that helps at the stall 
force is also able to make the completion of backward steps rather unlikely in the 
unloaded situation. 



4 Conclusions 



As the results indicate kinesin under no load could march along MT in the for- 
ward direction even without neck linker docking. Under the stall force, however, 
neck linker docking seems to be crucial to forward stepping, as for A'd = (cor- 
responding to no docking) the characteristic forward binding time tf is way over 
the desired value of 0.3 s. If the binding rate constant fcb near the MT is not ex- 
pected to be larger than its bulk value of 20 s^^fiM^^ by more than a factor of 



3, then we get very strict constraints for the neck linker: it should contain at least 
A^ = 6 Kuhn segments, and dock all along the head with a horizontal projection of 
Ld = 3.5 nm. Moreover, since the value of r quickly decreases as N is increased, 
the best strategy for kinesin to walk forward upto a loading force of —7 pN is to 
have a neck linker of about A^ = 6 Kuhn segments (corresponding to 13-14 amino 
acids), with N^t, = A or 5 segments capable of docking along the head. These 
values are very close to those known from the structure of kinesin, suggesting that 
kinesin is a very finely tuned motor protein whose parameters are set to live up 
to the highest load possible. This coincidence not only validates our calculations, 
but also justifies the relevance of neck linker docking in the motility of kinesin. 

5 Appendix 

Here we derive the formula for p%{R), the three dimensional probability density 
of the end-to-end distance of a FJC with A^ segments at zero loading force. Since 
N can be quite small we are not allowed to use the long FJC approximation of a 
polymer chain. 

Let p{R) be any three dimensional probability density with spherical sym- 
metry so that i? = is the origin. Then the relationship between p{R) and its 
projection p{x) to the x axis is: 

1 dp(x 



p{R) 



2n\R\ dx 



(11) 

R\ 



Since the end of a rigid rod can only move on the surface of a sphere the three 
dimensional density for an A^ = 1 FJC is a Dirac delta function on the surface 
of a sphere with radius Ik- Pi(-R) = ^(/k — |-R|)/(4/k^). By substituting p?(-R) 
into Eq. (11) and integrating it we get for Pi{x) a constant function with a value 
of 1/(2/k) and a support between —Ik and Ik- After autoconvoluting this function 
A^ times to get p%{x), and then converting it back to p%{R) with the help of Eq. 
(1 1), we arrive at the desired formula. 
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